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Abstract 

We show that the arguments forming the basis for the claim that the con- 
ditions for total neutrino conversion derived and studied in detail in jl], ^] "are 
just the conditions of the parametric resonance of neutrino oscillations supple- 
mented by the requirement that the parametric enhancement be complete", 
given in Q have flaws which make the claim physically questionable. We show 
also that in the case of the transitions in the Earth of the Earth-core-crossing so- 
lar and atmospheric neutrinos the peaks in the relevant transitions probabilities 
P a b, associated with the new conditions, max P a b = 1, are of physical relevance 
- in contrast to what is suggested in Q. Actually, the enhancement of P a b in 
any region of the corresponding parameter space are essentially determined by 
these absolute maxima of P^. We comment on few other aspects of the results 
derived in ||, [| which have been misunderstood and/or misinterpreted in 
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1. In a Comment on our results derived and discussed briefly in JT| and in detail in 
and on the results derived by one of us in ||, Akhmedov and Smirnov claim in 
particular, that "the conditions for total neutrino conversion studied by Chizhov and 
Petcov [|]]] Q are just the conditions of the parametric resonance of neutrino oscillations 
supplemented by the requirement that the parametric enhancement be complete." We 
show below that the arguments forming the basis for these claims have flaws which, 
in our opinion, make the claims physically questionable. We show also that in the 
case of the transitions in the Earth of the Earth-core-crossing solar and atmospheric 
neutrinos the peaks in the relevant transitions probabilities P Q/ g, associated with the 
new conditions, max P a p = 1, are of physical relevance - in contrast to what is 
suggested in M. Actually, the enhancement of P a p in any region of the corresponding 
parameter space is essentially determined by these absolute maxima of P a p. We 
comment on few other aspects of the results derived in ^ [2|, which have been 
misunderstood and/or misinterpreted in [|]]. 

The form of the new conditions for a total neutrino conversion in a medium con- 
sisting of two or three (nonperiodic) constant density layers, derived in [p], §, the 
region of the parameter space (i.e., the Am 2 /E — sin 2 29 plane) where they can be 
realized, and the physical interpretation of the corresponding absolute maxima of the 
neutrino transition probabilities of interest as being caused by a maximal constructive 
interference between the amplitudes of the neutrino transitions in the (two) different 
constant density layers, found in our studies, made us conclude in [[[], || that these 
new conditions differ from the conditions for parametric resonances in the neutrino 
transitions, discussed in the articles |], |7j and possible in a medium with density, 
varying periodically along the neutrino path. The Comment M does not provide 
viable arguments against our conclusions. 

2. Consider transitions of neutrinos crossing a medium consisting of i) two layers 
with different constant densities N\ >2 and widths Li,2, or of ii) three layers of constant 
density, with the first and the third layers having identical densities Ni and widths 
Li, which differ from those of the second layer, N 2 and L 2 0, 0, 0. Suppose the 
transitions are caused by two-neutrino mixing in vacuum with mixing angle 9. Let 
us denote by 9i and 20 i; i = 1,2, the mixing angle in matter in the layer with density 
Ni and the phase difference acquired by the two neutrino energy-eigenstates after 
neutrinos have crossed this layer. It proves convenient to introduce the quantities 
(see, e.g., @): 

cos $ = Y = c x c 2 - cos(29 2 - 29 1 )s 1 s 2 , (1) 
X 2 = 1 - F 2 , (2) 
n% sin $ = X 3 = — (s\c 2 cos 29\ + c\s 2 cos 2^2), (3) 

2 By using the verb "studied" in connection with our new conditions for a total neutrino conversion 
in a medium jl], ||] , the authors of Q suggest indirectly that these conditions were already considered 
in the literature before the publication of our papers. We would like to note that none of our new 
conditions for a total neutrino conversion in a medium were derived, postulated and/or discussed 
in some form in the literature on the subject published before the articles |l| |). 
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X3 being the third component of the vector [] X = (X\, X2, X3), whose first two 
components are also given in terms of 0j and fa (see 0). The probability of the 
transition v a — > Uf, (i.e., v e — > ^( T ), — > v e , etc.) after neutrinos have crossed n 
alternating layers with densities N\ and N2 is given according to @] by: 

P{y a - z/ b ; nL) = ( 1 - J sin 2 $ p = ^ + 2 ^ 2 sin 2 $ p , (4) 

where $ p = (n/2)$ if the number of layers n is even, and 

$ p = + f , <P = arcsin f Sl sin 26 1 /^l - X|/|X| 2 ^) (5) 

for odd number of layers with the first layer having density N]_. 

The first thing to note is that the expression for ip in eq. (|5|), given in ||, is 
strictly speaking incorrect: it is valid only if Z > 0, where 

Z = s 2 sin 2^ + si(sin20i)y. (6) 

The correct expression for (p for arbitrary sign Z reads: 

<p = arctan(sisin(20 1 ) \X\/Z) . (7) 

The authors of |4j] demonstrate the same imprecision in eq. (2) of their Comment 
the functions arccosF and arcsin |X| have different defining regions and it is incorrect 
to write arccos V = arcsin |X|: this equality is obviously wrong when Y < 0. 

According to the authors of "Eqs. @ and fl5|) describe the parametric oscil- 
lations with the pre-sine factor in (|]) and $ p being the oscillation depth and phase." 
and further "Parametric resonance occurs when the depth of the oscillations becomes 
equal to unity. The resonance condition is therefore X% = 0." 

In the two-layer case i) considered by us in , which is relevant for the present 
discussion, one has <3> p = $ (n = 2). This result and eqs. flU) and (@) imply that 
actually sin 2 $ p = Xf + X\ + X\ . Correspondingly, the "parametric-resonance" form 
in which the authors of 0] cast the probability P{y a — > u^; nL) in this case is artificial: 
the probability is given by 

P{u a -> v h - 2L) = X\ + Xj = l-Y 2 - Xl (8) 

where we have used eq. @. Therefore any resonance interpretation of the probability 
P{y a — > v\y\ 2L) based solely on the eq. (f|) with n = 2 seems to us physically 
unjustified. The new conditions for a total neutrino conversion in a medium follow 
in this case from the form (^) of P{y a — > v h \ 2L) and read JT|, JJ: 

Y = 0, X 3 = 0. (9) 

3 To facilitate the understanding of the main points of our criticism of g, we use the same 
notations as in 0] for most of the quantities discussed. 
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It should be clear from eqs. (§) and (^j that the condition X3 = alone does not 
ensure the existence even of a local maximum of the probability P(u a — *> u^; 2L). 

It is not guaranteed a priori that the equations in ([|) have non-trivial solutions 
in general, and in the specific case of transitions of neutrinos crossing the Earth core 
and the Earth mantle on the way to the detector. As we have shown in []1], 0, the 
solutions of the conditions (P) i) exist and are given by 



solution : 



tan 0i = ±^ 
tan 02 = ±1 



cos(2fl 2 



cos(20i) cos(2fl 2 - 20i) ' 

/ -cos(2fli) 
cos(2# 2 ) cos(2# 2 - 20i) ' 



where the signs are correlated, and ii) they can be realized only in the region deter- 
mined by the three inequalities: 

( cos 26 1 > 

region cos 2^2 < (11) 

[ cos(2# 2 - 20i) > 0. 

It was demonstrated in [l], [| as well that the two conditions in eq. (fj), or equivalently 
the solutions expressed by eqs. ( [T0| ) and (pi), are conditions for a maximal construc- 
tive interference between the amplitudes of the neutrino transitions in the two layers. 
Thus, a natural physical interpretation of the absolute maxima of P[y a — > z/&; 2L) 
associated with the conditions @ is that of constructive interference maxima. 

It should be clear from the above arguments that we do not see physical reasons 
to call X3 = in the case under discussion a "parametric resonance condition". 
Using the trick of |4[] one can easily cast the probability of two-neutrino oscillations 
in vacuum and in matter with constant density, for example, in the "parametric- 
resonance" form (HI). The analog of the condition X 3 = reduces in these cases 
respectively to sin 2 29 = 1 and sin 2 2# m = 1, 9 m being the mixing angle in matter. 
Thus, according to the terminology suggested in ||], one should call the condition 
sin 2 26 = 1 and the MSW resonance condition sin 2 29 m = 1 "parametric-resonance 
conditions". One can use such a terminology, of course, but this is not justified by 
the physics of the process and nobody uses it. The situation in the two-layer case 
discussed above is in essence the same. 

Analogous results and conclusions are valid in the case of medium with three 
layers (case ii)) considered in 0, [I], {n = 3 in eqs. (|) and (|)). Using the correct 
expression for <p, eq. (|7[), one finds again that the "parametric resonance" form in 
which the authors of 0j| write the probability of interest P{y a — > v b\ 3L) is artificial: 
the probability has the same form as in eq. (||) @ 

P(iy a ^u b ; 3L) = 1-Y 2 -Xl (12) 

where 

Y = -c 2 + 2 Cl F, (13) 

and 

X 3 = - S2 cos 2# 2 - 2s a cos(20i) Y. (14) 
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The conditions for a total neutrino conversion in this case read 

Y = 0, X 3 = 0. 



(15) 



Their solutions in the case of the inequality iV\ < N 2 (corresponding to the relation 
between the densities in the Earth mantle and core) were given in |T|, |2| and have the 
form: 



solution : 



tan <pi 
tan <p2 



' - cos 26 2 
cos(2# 2 -40i)' 
cos 2#i 



(16) 



cos(2# 2 ) cos(2fl 2 -40i) 
where the signs are again correlated. The solutions can only be realized in the region 



region 



A®: 



cos(2fl 2 ) < 0, 
cos(2fl 2 - 40i) > 0. 



(17) 



It is easy to show that for any number of layers n, the denominator in (|4]) is 
always canceled by sin 2 <3> p and P(v a — > v h \ nV) is just a polynomial without any 
resonance-like feature. Indeed, for even n, as can be shown, we have 



sin$ p = sin$ U n / 2 -i(cos $), 

where U n (x) is well-known Chebyshev's polynomial of the second kind 
of odd number of layers (n > 3) one finds 



P(v n 



u b ; nL) 



si sin 29 1 cos ^ - - + ZU n -3 (cos $) 



In the case 
(19) 

Finally, similar considerations apply to the probability of the v 2 — > z/ e transitions, 
^2 being the heavier of the two vacuum mass-eigenstate neutrinos, in the three-layer 
case ii), P{y<i — > 3L). This probability can be used to account for the the Earth 
matter effects in the transitions of solar neutrinos traversing the Earth: P{yi — > 
v e \ 3L) corresponds 0] to the case of solar neutrinos crossing the Earth mantle, the 
core and the mantle again on the way to the detector. As was shown in ||, the 
conditions for a total u 2 — > v e conversion, max P{y 2 — > v e ] 3L) = 1, read: 



Y = 0, X' = 0, 



where Y is given by eq. (p!3|) and 

= -s 2 cos(20 2 



2si cos(20! 



(20) 



(21) 



The solutions of the conditions ( |2"0"D providing the absolute maxima of P{v 2 v e \ 3L) 
and the region where these solutions can take place were given in |l], 0; they can 
formally be obtained from eqs. (IB) and ( |T7D by replacing 26*! and 29 2 with (29 1 — 9) 
and (202 - 9). 

The three sets of two conditions, eqs. ©, flTBD and (^0| ), and/or their solutions 
(e.g., eqs. flip]) and (0)), and/or the regions where the solutions can be realized (e.g., 
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eqs. ( |TI] ) and (JT7\)), were not derived and/or discussed in any form in H ||, [7j or in 
any other article on the subject of neutrino transitions in a medium published before 
0, None of them follows from the conditions of enhancement of P{y a —>■ ^b) found 
in H ||, ^] and thus they not a particular case of the latter. That is the reason we 
used the term "new conditions for a total neutrino conversion in a medium" for them. 

The authors of |4| write further: "One well known realization of the parametric 
resonance condition", i.e., of X 3 = 0, is |, 0, [| f] "ci = c 2 = 0, or 

20! = 7T + 27th, 20 2 = vr + 2nk 2 , fci, fc 2 = 0, 1, 2, (22) 

independently of the mixing angles." Contrary to what the authors of claim, the 



two conditions in eq. (|22| ) were not given in the articles || g, 0: what one finds in 
these articles at most is the condition 20i + 202 = 2tt + 27rA; which is not equivalent 
to the two conditions in eq. (p2[). The two conditions in eq. ( p2|) were discussed in 
detail for the 3-layer case in j3|. Moreover, as we have shown in || and would like 
to emphasize here again, the conditions c\ = 0, c 2 = fry themself do not lead to 
a maximum of the neutrino transition probabilities of interest in the neutrino energy 
variable, unless a third nontrivial condition is fulfilled. This third condition has the 
following form for the v a — > (i.e., v e — > ^ M ( T ), — *> z/ e , etc.) transitions in the 
two-layer (n = 2) and three- layer (n=3) medium cases i) and ii), respectively [[|: 
cos(2# 2 - 2Q X ) = and cos(2# 2 - 40i) = 0. For the probability P{v 2 u e ; 3L) 
it reads ||: cos(2# 2 — 46i + 6) =0. It is not difficult to convince oneself that the 
indicated sets of three conditions represent possible solutions respectively of ([|), ( |T5D 
and (^) More generally, the condition X3 = alone does not guarantee the 
existence even of a local maximum of the neutrino transition probabilities of interest. 

In what regards the article by Q.Y. Liu and A. Yu. Smirnov quoted in ||] in con- 
nection with the conditions (p2|) (see ref. [5] in 0]), these authors noticed that in the 
case of muon neutrinos crossing the Earth along the specific trajectory characterized 
by a Nadir angle h S 28.4°, and for sin 2 26 S 1 and Am 2 /E = (1-2) x 10" 4 eV 2 /GeV, 
the —* v s transition probability, v s being a sterile neutrino, is enhanced. The au- 
thors interpreted this enhancement as being due to the conditions 2<pj = it, j = 1, 2, 
which they claimed to be approximately satisfied. Actually, for the values of the 
parameters of the examples chosen by Q.Y. Liu and A. Yu. Smirnov to illustrate 
their conclusion one has 2(f>i = (0.6 — 0.9)7r and 2<p 2 = (1.2 — 1.5)7r. The indicated 
enhancement is due to the existence of a nearby total neutrino conversion point 
which for h = 28.4° is located at sin 2 26 S 0.94 and Am 2 /E = 2.4 x 10" 4 eV 2 /GeV 
and at which 2<pi = 0.9tt and 2<p 2 = Lin. We have also found in 0] that for each 
given h ^ 30° there are several total neutrino conversion points at large values of 
sin 2 26, at which the phases 20i and 202 are not necessarily equal to n or to odd 
multiples of it (see Table 3 in 0). Thus, the explanation of the enhancement offered 
by the indicated authors is at best qualitative and incorrect in essence. 

That X 3 = alone is a condition for local maxima of the probability P[y a — > 
u b ; nL) was suggested in E. Kh. Akhmedov, Nucl. Phys. B538 (1999) 25 (to 
be quoted further as NP B538, 25), on the basis of the form of the probability in 

4 We quote here only the references which, in our opinion, are relevant for the present discussion. 
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eq. (ffl). However, as we have already emphasized, the condition X3 = alone 
does not guarantee the existence even of a local maximum of the neutrino transition 
probabilities P{y a — > Uf,; 2L) and P[y a — > v\,\ 3L) . This should be clear from the 
"natural" expressions for the probabilities P{y a — > IV) and P{y a — > z/j,; 3L), 
given by eqs. (H) and (^). Of the three solutions for the extrema of P{y a — > v\,\ 3L) 
found in NP B538, 25, the solution c\ = 0, c 2 = was already discussed in detail in 
ref. m (compare eqs. (11) - (16) and (24) in || with conditions (1) on page 37 of 
NP B538, 25), while the other two correspond to MSW transitions Q (they were also 
briefly discussed in ||). 

4. The authors of claim that "The existence of strong enhancement peaks 
in transition probability P rather than the condition P=l is of physical relevance.", 
although they do not give an example of a relevant strong enhancement peak (i.e., 
local maximum). As we have shown in |2| (see also [[|), the solutions given by eqs. (10) 
and flTB] ) and those for the probability P{y 2 v e] 3L) are realized in the transitions 
in the Earth of the Earth-core-crossing neutrinos (solar, atmospheric, accelerator) and 
lead to observable effects in these transitions. From the extensive numerical studies 
we have performed in the realistic cases of transitions of neutrinos in the Earth (e.g., 
neutrinos crossing the Earth core on the way to the detector) we do not have any 
evidence about the presence of significant local maxima in the neutrino transition 
probabilities of interest not related to the peaks of total neutrino conversion. Actually, 
our studies show that only the peaks of total neutrino conversion are dominating in 
P(y a ~^ v b\ 2L), P(u a — > u b ; 3L) and P{y2 —* v e \ 3L), and correspondingly determine 
the regions where these probabilities can be significant in the corresponding space of 
parameters. The peaks considered, e.g., in |3j (see figs. 1 and 2) are points on the 
"ridges" formed by local maxima, e.g., in the energy variable at fixed values of the 
other parameters, leading to the peaks of total neutrino conversion, discovered by us 
(see figs. 6-9 and Tables 5 - 6 in f2|). As one can convince oneself using Figs. 6-9 
from [Q], all maxima in Figs. 4, 5, 10 - 13 and 15 in ||, including the relatively small 
local ones, are related to (and determined by) the presence of corresponding points 
(peaks) of total neutrino conversion. 

5. The authors of 0, [7) studied the effects of small density perturbations on the 
neutrino oscillations, while we have investigated in |3], III, 0] the different physical 
problem of large "perturbations" of density. In || neutrino oscillations in a medium 
consisting of even number n of alternating layers with densities Ni and N 2 have been 
considered. However, i) the two layers were assumed to have equal widths, Li = L 2 , 
and i) an enhancement of the neutrino transitions was found to take place for small 
vacuum mixing angles at densities Ni^ much smaller than the MSW resonance den- 
sity. The authors of H |6| were interested in and found the conditions for the classical 
parametric resonance in neutrino oscillations, which can take place after many periods 
of density modulation of the oscillations. From reading the articles || |6], 0] it becomes 

5 Let us note that in what regards the cases of neutrino transitions studied in [|[ |l| ||], the article 
NP B538, 25, contains a rather large number of incorrect statements and conclusions. Most of these 
statements are concentrated in Section 4 of NP B538, 25, where the author discusses the realistic 
case of transitions of neutrinos crossing the Earth core we were primarily interested in in [l]. |^] . 
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clear that their authors had in mind astrophysical applications of their results [], and 
the results obtained in [§, ^ |7| may indeed have astrophysical applications. Our stud- 
ies |3|, [I], were concerned primarily with the neutrino oscillations in the Earth. As 
a consequence, the conditions of the enhancement of P{y a — > v^) obtained in || |], [7| 
differ from those found by us in |3|, [I], In the two- and three-layer medium cases i) 
and ii) discussed by us (e.g., in the case of the Earth) the enhancement found in ||, 
for example, is not realized. 

To summarize, after studying ref. Q and the references quoted therein one can 
conclude that the new conditions for a total neutrino conversion in a medium found 
and studied in []!], (i.e., the three sets of two conditions, eqs. ([|), (|1^) and fl20D , 
and/or their solutions (e.g., eqs. (|T0D and (|16D ), and/or and the regions where the 
solutions can be realized, eqs. ([11]) and (0)), were indeed new: they were not derived 
and/or discussed in any form in |5|, ||, 0] or in any other article on the subject of 
neutrino transitions in a medium published before |]I|, |2[. None of them follows from 
the conditions of enhancement of P{y a — > v b ) obtained in || |6|, [7j and thus they are not 
a particular case of the latter. In M , in particular, the derivation of these conditions, 
first given in ||, is reproduced. Most importantly, the new conditions for a total 
neutrino conversion were shown in || (see also |J) to be realized in the transitions in 
the Earth of the Earth-core-crossing neutrinos (solar, atmospheric, accelerator) and to 
lead to observable effects in these transitions - contrary to the claims made in Q. As 
for the physical interpretation of the associated new effect of total neutrino conversion 
in the cases of two-layer and three-layer medium we have considered, we have proven 
that this is a maximal constructive interference effect. The interpretation of the effect 
based on the expression ([|) for the neutrino transitions probabilities, offered in 0, as 
we have pointed out, is not convincing: expression (4), in particular, has an artificial 
resonance-like form. For the studies of the new effect of total neutrino conversion in 
the Earth, eqs. (H) and ( |i~2|) represent one of the several possible natural expressions 
for the relevant neutrino transition probabilities. The rest is terminology. 



6 In all the concrete examples and corresponding plots in Q, for instance, the physical quantities 
having a dimension of length are given in units of the solar radius and are comparable with it; the 
values of the densities used in the examples differ substantially from those met in the Earth - they 
are noticeably smaller than the Earth mantle and core densities. 
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